Lagrangian statistics in forced two-dimensional turbulence 
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We report on simulations of two-dimensional turbulence in the inverse energy cascade regime. 
Focusing on the statistics of Lagrangian tracer particles, scaling behavior of the probability density 
functions of velocity fluctuations is investigated. The results are compared to the three-dimensional 
case. In particular an analysis in terms of compensated cumulants reveals the transition from a 
strong non-Gaussian behavior with large tails to Gaussianity. The reported computation of correla- 
tion functions for the acceleration components sheds light on the underlying dynamics of the tracer 
particles. 
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Introduction In recent years, the Lagrangian descrip- 
tion of turbulent flows has attracted much interest from 
the experimental point of view [l|, Q as well as in numer- 
ical [1^ Q and analytical investigations Q. This is not 
only due to the relevance of the Lagrangian approach 
for applications like turbulent mixing and the dispersion 
of pollutants. But fundamental turbulence research also 
benefits from this alternative description and its rela- 
tion to the Eulerian formulation. The turbulent state, 
for example, is completely represented by the accelera- 
tion of the tracer particles. The classical Kolmogorov 
theory (K41) predicts self-similar scaling of the proba- 
bility density functions (pdfs) of the velocity increments 
in homogeneous isotropic turbulent flows. It is a well 
known fact that in three dimensions the assumption of 
self-similarity is violated for Eulerian velocity increments. 
This is referred to as Eulerian intermittency. We know 
from recent experiments and numerical simulations, that 
intermittency is observed in the Lagrangian picture as 
well. For two-dimensional turbulence the situation is 
different: although numerical investigations show devi- 
ations from Gaussianity Q for the Eulerian velocity in- 
crements, the inverse cascade provides an Eulerian flow 
field with scaling properties compatible with the K41 pre- 
dictions. For this reason, the inverse energy cascade in 
two-dimensional turbulence seems to be the ideal system 
to shed light on the Lagrangian statistics and its relation 
to the Eulerian case. 

The paper in hand presents results from numerical in- 
vestigations of the Lagrangian dynamics of forced two- 
dimensional turbulence. Thereby we focus on intermit- 
tency. The remainder of this article is structured as fol- 
lows. After briefly discussing the system under considera- 
tion and summarizing some facts regarding the numerics, 
we present results on the Eulerian increment statistics. 
The main part of the article is devoted to the study of 
the Lagrangian increment pdfs. Finally we discuss prop- 
erties of the acceleration correlation. 

Basic equation The behavior of two-dimensional fluid 
motion is governed by the vortex transport equation 



dfUJ = —u ■ Voj + v/S.Lj — juj — A/, 



(1) 



= —dxx4> are connected to the vorticity field u via 
the stream function cf) by = —lo. We investigate 
two different types of forcing /. One with a rapidly 
decaying spatial correlation function {f{x + r)f{x)) ^ 
exp(— r^/2Z^) [5] where Ic is the length-scale of the cor- 
relation. The second is confined to a small shell of 
wavenumbers in Fourier space 0] leading to a long spa- 
tial correlation. Both forcings are 5-correlated in time 
and have the property to inject energy at small scales 
into the system. The damping term — 7W in ([T]) extracts 
energy at large scales from the system and avoids the 
generation of a large scale flow. 

Numerical method The integration of ([T]) was performed 
by a fully dealiased pseudo-spectral method on a doubly 
periodic square domain with side length 27r and 1024^ 
grid points. For numerical reasons, the viscous term 
is replaced by a hyperviscous term of order eight. The 
time evolution for the field and the tracer particles was 
achieved by a fourth order Runge-Kutta scheme. Af- 
ter reaching stationarity we introduced 10^ particles with 
random but equally distributed initial positions into the 
flow and monitored their trajectories for about 200 La- 
grangian integral times T/. The velocity and the acceler- 
ation of the tracers were determined by a bicubic interpo- 
lation scheme. Except for Fig[T] and Fig[7] all presented 
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with the vorticity lu = uj{x,t) and the velocity u — 
u{x,t). The velocity components u^-^ = dx24' and 



FIG. 1: Energy spectra E(k) of simulations using the forcing 
with the short (upper curve) and the long (lower curve) cor- 
relation length. The line denotes the Kolmogorov prediction 
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FIG. 2: Rescaled pdfs of the Eulerian velocity increments for 
the distances r — 0.08,0.1,0.2,0.6. For comparison a Gaus- 
sian pdf is shown. Small deviations from Gaussianity exist 
and can be quantified by the moments 0]. 
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FIG. 3: Rescaled pdfs of Lagrangian velocity increments for 
the time lags r = 0.045, 0.09, 0.22, 1.79T/ (from outer to inner 
curves). The most inner curve is a Gaussian pdf. The pdfs 
are vertically shifted. 



results are obtained using the forcing with the rapidly 
decaying spatial correlation. 

Eulerian statistics In order to check the parameters for 
the numerical integration, we first analyze the Eule- 
rian velocity field. FiglT] shows the energy spectrum for 
both kinds of forcing together with a line showing the 
K41-scaling For a stationary velocity field, the 

longitudinal Eulerian velocity increments are defined as 
6ve{x, r) = Sve{x, r) f with Sve{x, r) ~ u{x + r) —u{x) 
and f = r/r. If we additionally assume isotropy and ho- 
mogeneity of the flow, we can write Sve{x,r) — Sveir). 
The pdfs p{Sve{r)) are scaled to unit standart deviation 
by ap{6ve{r)/a) with a = (5ue(r))^/^. FigE] shows the 
rescaled pdfs for different r. The shape does not vary 
with r and hence the pdfs are self-similar. This is in 
agreement with experimental Q and numerical Q stud- 
ies and leads to the conclusion that intermittency is ab- 
sent in the inverse energy cascade as far as the Eulerian 
increments are concerned. 

Lagrangian velocity statistics In the Lagrangian frame of 
reference the velocities are recorded along the trajecto- 
ries of tracer particles v{y, t) — [u{x, i)] j.^x(j/ 1) ' where 
y is the starting position of the tracer and X{y,t) is 
its current position. Velocity fluctuations are character- 
ized by the pdfs of the Lagrangian velocity increments 
Svi (r) = Vi{t + T) — Vi (t) with i = xi , X2 . Due to isotropy 
we average with respect to both spatial components. 
Therefore in the following the Lagrangian velocity incre- 
ment will be denoted as Sv{t). The moments of p{Sv{t)) 
are known as the structure functions 5'p(r) = {Sv{t)p). 
FigEl shows pdfs of the Lagrangian velocity increments 
for several time lags r. For time lags of the order of T/, 
the pdfs are close to a Gaussian distribution whereas for 
small r the pdfs show large tails. A rescaling of the dis- 
tribution functions resulting in a collapse to a universal 
distribution is not possible. Accordingly, classical Kol- 
mogorov scaling can not be observed in the Lagrangian 
frame in contrast to the Eulerian case. Deviations from 
the Gaussian shape can be quantifled by the compen- 
sated cumulants A„ = c„/(t". The c„ are the cumulants 



connected to the characteristic function C{k) of a pdf by 



C{k) = exp 



E'--n 
7^ 



(2) 



and (7^ is the variance (corresponding to C2). 

For a Gaussian distribution all A„ with an order 11 
higher than 2 vanish. The compensated cumulants can 
easily be computed from the structure functions Sp. For 
symmetric pdfs, A4 = S4/S2 — 3 is the excess kur- 
tosis and the sixth order normalized cumulant reads 
As = (^6 - lhSiS2)/Sl + 30. 

In FigH] we see a log- log plot of A4 and Ag . The crucial 
point is that for a self-similar signal the kurtosis should 
be a constant at least in the region where self-similarity 
of the pdfs is expected. In our simulations the kurtosis 
depends strongly on r and for intermediate times follows 
a power-law as can be seen in FigHl This also holds for 
Ag which however decays faster than A4. As a reference 
we plotted two lines following the power laws r~^'^^ and 
T^^-^. During the decay of the cumulants the pdfs con- 
verge to the Gaussian shape. In order to investigate the 
universality of the observed behavior we used the data 
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FIG. 4: Compensated cumulants of order four (lower curve) 
and six (upper curve) for 2D turbulence. As a guide for the 
eye t~^'^^ and r~^'® are shown. 
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FIG. 5: Compensated cumulants of order four (lower curve) 
and six (upper curve) for 3D turbulence together with r~^'^^ 



and r 
time r„ 



The time is given in multiples of the Kolmogorov 



FIG. 6: ESS plot of (lower curve) and 5*5 (upper curve). 
The lines show the ESS scaling laws with exponents (lines 
ranging from 5| = 0.04 to Sa = 0.4) and Cp (from S2 = 0.015 
to S2 = 0.9) obtained by the two fitting procedures (see text). 
The dashed lines correspond to the K41 scaling. 



provided by to calculate the same quantities for 

three-dimensional turbulence (see Fig[5]). Again scaling 
behavior can be detected for intermediate times and for 
comparison we added power-laws with the same expo- 
nents as in Fig[3) For very small and very large r the 
shape of the compensated cumulants show differences be- 
tween 2D and 3D. We devote this to the fact that in 2D 
the energy is injected on the small scales and mainly dis- 
sipated at the large scales whereas in 3D the situation is 
complementary. Again our results give strong evidence 
for intermittency in two-dimensional Lagrangian turbu- 
lence. 



(Figl?]) as well as for the cumulants are qualitatively and 
quantitatively similar to the situation with short spatial 
correlation. This leads to the conclusion that the ob- 
served deviation of the Lagrangian pdfs from the K41 
prediction is very robust and seems to be independent of 
the type of forcing. 

Acceleration correlations The path of a Lagrangian tracer 
particle starting at the position y is uniquely defined by 
the acceleration acting on the particle. The acceleration 
is given by the right-hand side of the Navier-Stokes equa- 
tion 
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0.557 ± 0.002 


1.267 ± 0.007 


1.35 ± 0.018 


1.313 ± 0.033 
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0.557 ± 0.003 


1.313 ± 0.008 


1.45 ± 0.019 


1.588 ± 0.029 



TABLE I: ESS scaling exponents for 2D turbulence 



Scaling of the Lagrangian structure functions Addi- 
tionally to the cumulants, the structure functions are 
computed to characterize the pdfs. As no scaling re- 
gion is visible for the structure functions we have to rely 
on the Extended Self Similarity (ESS) technique [10| in 
order to estimate scaling exponents. To apply ESS we 
have to use the structure functions for the absolute values 
of the velocity increments S*{t) = {\Sv{t)\p). Standard 
arguments of dimensional analysis lead to the scaling be- 
havior S* ~ r^p with (p = p/2. The ESS-plot is shown in 
FiglG] Estimating the exponents up to order five in the 
spirit of [1] leads to the values shown in TablH We 
also performed the analysis for larger values of 5*2 (corre- 
sponding to larger t) 0] resulting in the exponents Cp- In 
both cases the exponents deviate strongly from the K41 
predictions which is in agreement with the observation 
that the excess kurtosis is not constant. 

Dependence on the forcing To study the effect of the dif- 
ferent forcings on Lagrangian observables, we also per- 
formed simulations with a forcing limited to a small num- 
ber of Fourier modes. In this case, the results for the pdfs 



a{y,t) 

= [~Vp{x, t) + v^u{x, t) + V X f{x, . 

(3) 

where the pressure p{x, t) is related to the vorticity by 
Ap = V • [ti X w] — ^Aii^. Besides this fact, the accel- 
eration is also of central interest for turbulence model- 
ing prl. [T^. Particularly in two dimensions, it is conve- 
nient to split up a(x,t) into a component parallel and 
a component perpendicular to the current velocity of 
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FIG. 7: Rescaled pdfs of Lagrangian velocity increments for 
the same time lags as in Fig|3] In this case the forcing was 
confined to a small number of Fourier modes. 
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FIG. 8: Correlation functions Ca^ij) (upper curve), Ca||(r) 
(lower curve) and C-vij) (middle curve). The inset presents a 
log-log plot of —Cay (r). The line corresponds to t~^'^ . 

the tracer. The latter is sensitive to circular motions. 
The two components are defined as a\\ = a ■ u/\u\ and 
a± = a ■ {—Uy,Ux)/\u\. As pointed out in long cor- 
relation times of the velocity increments and the accel- 
eration play a key role for the occurrence of Lagrangian 
intermittency. This poses the question whether the ob- 
served deviation from self-similarity of the increment 
pdfs in 2D is also connected to long correlation times 
of the acceleration. In FiglS] we can see the correlation 
functions C/i(t) with h — aj_,a||,i> for the two compo- 
nents of the acceleration and for the velocity defined as 
Chir) = {h{t)h{t -h r))/(/i2(t)). The correlation for the 
parallel component decays very fast and approaches zero 
after passing a minimum at negative values. The per- 
pendicular acceleration component also decorrelates very 
fast for small r. For r bigger than the time correspond- 
ing to the minimum of Cay (t) it bends off into a region 
where it decays much slower leading to a very long cor- 
relation time. The same behavior is also observed for 3D 



turbulence [l|, [3, [TB| . In [TJ] it was related to the spi- 
raling motion in a vortex filament. Here we also observe 
events where the particle runs through loops. Because 
the motion is confined to a plane, there is no movement 
in the third spacial dimension that could contribute to 
the decorrelation of a±_ . For comparison also C„ is shown 
in Figl51 The inset in the same figure demonstrates that 
Cay follows a power-law in the range between its mini- 
mum and Tf. The results for the temporal correlations 
of the acceleration suggest that the stochastic process for 
the velocity increments is essentially non-Markovian, as 
has been emphasized in [l 6| . 

Conclusion We presented a detailed investigation regard- 
ing the statistics of tracer particles in the inverse energy 
cascade regime of two-dimensional turbulence. For differ- 
ent types of forcing, we detect the same deviations from 
self-similarity for the Lagrangian velocity increment pdfs. 
This is strong evidence in favor of Lagrangian intermit- 
tency in the inverse energy cascade. It is of particular 
interest as for the Eulerian frame no intermittency can 
be detected. Any attempt to relate the two frames of ref- 
erence has to incorporate this fact. The observation that 
in 2D and 3D the compensated cumulants show the same 
scaling behavior for intermediate time lags suggest that 
the underlying dynamical process exhibits a certain de- 
gree of universality independent of the dimension. This 
view is supported by the fact that the acceleration com- 
ponents show long time correlations which are similar to 
the 3D case. The explanation for the scaling of the com- 
pensated cumulants remains an open question. 
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